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Abstract 

This paper investigates the impact of control field noise on the optimal manipulation of quantum 
dynamics. Simulations are performed on several multilevel quantum systems with the goal of 
population transfer in the presence of significant control noise. The noise enters as run-to-run 
variations in the control amplitude and phase with the observation being an ensemble average over 
many runs as is commonly done in the laboratory. A genetic algorithm with an improved elitism 
operator is used to find the optimal field that either fights against or cooperates with control field 
noise. When seeking a high control yield it is possible to find fields that successfully fight with the 
noise while attaining good quality stable results. When seeking modest control yields, fields can 
be found which are optimally shaped to cooperate with the noise and thereby drive the dynamics 
more efficiently. In general, noise reduces the coherence of the dynamics, but the results indicate 
that population transfer objectives can be met by appropriately either fighting or cooperating with 
noise, even when it is intense. 
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I. INTRODUCTION 



Control of quantum processes is actively being pursued theoretically 



experimentally 



3, 



DD 



and 



In practice control field noise and environmental interactions inevitably 
will be present, and the general expectation is that their involvement will be deleterious to- 
wards achieving control [5]. The present paper only considers the influence of field noise 
upon the controlled dynamics with the noise described by shot-to-shot pulse variations, as 
in the signal average experiments. A separate work will explore analogous issue for environ- 
mental interactions. Recent studies have considered several aspects of the influence of laser 



noise 



SB 



id, and this work aims to further explore the issue. Operating under 



closed loop|l2| in the laboratory will naturally deal with the noise as best as possible. Using 
simple models and simulations, this paper will explore the nature of robustness to noise 
and the possibility of cooperating with the noise rather than fighting against it, especially 
in the strong noise limit. Such a prospect has foundation in analogous classical stochastic 
resonance phenomena [l^] . 

The presence of field noise in quantum control is generally viewed as problematic, 
but recent successfu 



tolerated jl 
dynamics [20, 



15, 
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21. 
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experiments at least support the point that some noise can be 
jjl Irsl Irgj ] and that it may have a constructive effect in quantum 



]. A recent theoretical analysis |7J] on the impact of field noise upon 
optimal control indicated that an inherent degree of robustness can be anticipated by virtue 
of controlled observable expectation values being bilinear in the evolution operator and its 
adjoint. Some simulations of closed loop experiments also show that robust control is pos- 



sible, and that properly designed control fields can fight a gain st noise 



special circumstances, field noise even may be helpful 
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24 I25II . In some 



3 If ] . The pres- 



ence of field noise may also improve the convergence rate of the laboratory learning control 



processes |2j 

In this paper, we will investigate four model illustrations of optimal control in the presence 
of field noise, especially for cases of strong noise. It is naturally found that field noise is 
deleterious to achieving control if a high yield is desired, however we also find robust high 
yield solutions that are extremely stable with respect to noise. If a low yield is acceptable, 
then it is shown that field noise may even be helpful and the control field can cooperate 
with it. The present work primarily aims to demonstrate these findings in several model 
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illustrations. A follow on paper will analyze the mathematical and physical origin of the 
observed behavior j^. 

Section II will present the model and control formulation used in these studies. Section 
III describes the genetic algorithm employed to search for an optimal control field in the 
presence of noise, and some special algorithmic modifications are made for this purpose. 
Section IV presents a set of control simulations under various system and noise conditions. 
Some brief concluding remarks are given in section V. 



II. THE MODEL SYSTEM 



The effect of field noise on controlled quantum dynamics will be explored in the context 
of population transfer in multilevel systems characterized by the Hamiltonian H, 

H = H - uE(t) (la) 
H = Y,e v \v)(v\ (lb) 

V 

where \v) is an eigenstate of H and e v is the associated field-free eigen-energy. The noise-free 
control field E(t) is taken to have the form 

E (t) = S(t) A * cos M + °i) ' ( 2a ) 
1 

S(t) = exp [- (t - T/2) 2 /2a 2 ] , (2b) 

where {ui} are the allowed resonant transition frequencies of the system. The controls are 
the amplitudes {A®} and phases {0®}. 

In each simulation, a reference noise-free optimal control calculation is first performed 
with the cost function: 



J [E Q (t)] = \0 [E (t)} - T \ 2 + aF (3a) 

fo = EKf, (3b) 
l 

where Ot is the target value (expressed as a percent yield in the applications later) and 
O [-Eo(^)] is the outcome produced by the field E (t), and F is the fluence of the control field 
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whose contribution is weighted by the constant, a > 0. In the present work, O = \^f) {9f\ 
is a projection operator for the population of the target state |^/). 

Noise in the laboratory could take on various forms and arise from a number of sources^]. 
In keeping with laboratory practice, the achieved control will be measured as an ensemble 
average over the outcome of many control fields. Here the noise is modeled as run-to-run 
uncertainties in the amplitudes and phases in Eq. This work will take an extreme 

conservative view to explore the nature of achievable control when such noise is large and 
independent of each other in the different runs. The noise is simulated by introducing the 
parameters 7^ and 70, which are chosen randomly over the range [— 1^,1^] and [— Tg,Tg], 
respectively, 

4° = 4° + 72, <P = Q + '$- (4) 
What is referred to as "amplitude noise" alone actually allows for phase changes over the 
range [— r^lV]. By exploring the case of Tg 7^ 0, Ta = 0, the phase contribution will 
become clear. Different random values are selected for the amplitudes and phases over their 
respective ranges in replicate simulations i = 1, 2, ■ • ■ , M corresponding to an ensemble of 
M control experiments. The z'-th experiment is driven by field 

(t) =s{t)J2 A w cos Lit + e?A . (5) 

I 

and the net outcome of all the control experiments is the average, 

M 

:\ [ 



M 

<o (*)]>„ = ^EKV'([^(*)] J r)|^> 1 2 , (6) 



i=i 



where \ip ,T)) is the final state driven with the z-th field in Eq.fl^J). The bracket 

notation (£)] indicates a functional dependence on the field over the internal < t < T. 
The objective function in the presence of noise is: 

J N (A?, 0°) = |<0 -0 T \ 2 + aF . (7) 

In order to find a robust control field, for which the outcome is optimally insensitive to 
the noise, we incorporate the measurement standard deviation a into the cost function to 
finally produce: 

j(Al9f) = J N (l+(3a N [E^(t)]), (8) 
4 



where 

a N [E®(t)] = ^(0'[E^(t)}) N -(0[E(Ht)]) 2 N (9) 

and (0 2 ) N is defined similar to Eq.(jB|). The coefficient (3 > balances the magnitudes of 
the terms. Since J is minimized, pulses that produce a small variance in the objective are 
favored. 



III. A GENETIC ALGORITHM WITH AN IMPROVED ELITISM OPERATOR 
IN THE PRESENCE OF CONTROL NOISE 



Current quantum control experiments [12j employ genetic (GA) 



algorithms 



or evolutionary 



24, 



that GAs 



14j to identify optimal control fields. Simulations have shown 
can tolerate a high level of noise because they don't make sharp decisions, allowing for de- 
cision errors to be averaged out over the generations of the learning process. GAs can be 
very effective, but they also can be expensive to employ, especially when optimizing a shot- 
to-shot stochastic process because a large number of experiments are needed to perform the 
ensemble average. Many of the simulations in this work include significant field noise, and 
it was found that at least M = 10 4 runs are needed to get accurate ensemble averages. In 
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Here we 



order to speed up the simulations, many strategies have been proposed 
used a simple approach: only perform full statistical averaging with the most effective field 
at each generation, while for intermediate assessments of all other individuals use a small 
number of ensemble average runs. The detailed procedure to select an optimal field consists 
of the following steps: 

1. Initialize the GA. We assume there are two special individuals, one is denoted as 
" best" , which is the best member of population; another is denoted as " elitist" , which will 
be updated in each generation, elitist is the same as best when the GA starts. 

2. Perform the standard operations of reproduction, crossover and mutation to generate 
a new population for the next generation. Evaluate the fitness of each individual in the 
population with a small size ensemble (M < 10 3 in this work). 

3. Compare the objective values arising from elitist and best (the best member of the 
new generation): 

if J (elitist) < J (best) then re-evaluate the cost function of elitist with a large ensemble 
(M > 10 4 in this work); 



if J (elitist) > J [best) then re-evaluate the cost function of best with a large ensemble 
(M > 10 4 in this work); 

4. Compare the updated cost function value of elitist and best from the last step; 
if J (elitist) < J (best) then replace the worst member of population with elitist 

if J (elitist) > J (best) then set elitist = best; 

5. Go to step 2, and repeat the GA operations until convergence is fulfilled. 
Although there is always a chance that a good field is overlooked by this procedure, it 

proved to be effective in finding robust fields. The procedure is very efficient as it only 
performs extensive statistical averaging on the likely most effective fields. 

IV. NUMERICAL SIMULATIONS 

To demonstrate how quantum optimal control can either fight or cooperate with field 
noise, we performed four simulations with simple model systems. The first three simulations 
use the single path system in Figure 1(a) while the last simulation uses the double-path 
system in Figure 1(b). Equation (jHJ) is the cost function used to the guide GA choice of 
the control field. The average yields shown in all the tables typically involve an ensemble 
of M = 10 4 samples. The simulations used the Toolkit technique 37} to efficiently solve the 
Schrodinger equation. 

A. Model 1 

This model used the single path 5-level system in Figure 1(a) with eigenstates \i), i — 
0, • • • , 4 of the field free Hamiltonian Hq, having only nearest neighbor transitions with the 
frequencies Wqi = 1-511, U\ 2 = 1-181, uj 2 3 = 0.761, CU34 = 0.553 in rad-fs -1 , and associated 
transition dipole moments /x i = 0.5855, /ii 2 = 0.7079, // 2 3 = 0.8352, /i 34 = 0.9281 in 
1.0 x 10~ 30 C-m. The target time is T = 200 fs, and the pulse width in Eq.fjSJ) is s = 30 fs. 
The control objective is to transfer the population from the initially prepared ground state 
|0) to the final state |4), such that O = |4) (4|. As a reference control case, the algorithm 
is run first without any field noise using the cost function in Eq. (|3ajl . Figure 2 shows the 
optimal field for the two separate goals Ot = 100% and 2.25%. These two extreme cases of 
high and low target yield are chosen to illustrate their distinct control behavior in the present 
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of noise. The choice of the low yield target of 2.25% is arbitrary with similar behavior found 
for yields up to ~ 10%. In some applications even making a modest yield of a particular 
product or state for study would be significant. The physical criteria in the low yield limit 
is that sufficient noise be present to produce at least a small population in the target state. 
At the target time, the noise free fields in Figures 2(al) and 2(bl) drive 98.7% and 2.2% of 
the population to target state, respectively. 

In order to calibrate the strength of the noise, Table I shows the state population distri- 
butions when the system is only driven by noise of various amplitudes. When Ta > 0.07, 
the influence of the field noise is very large as it drives more than 40% population out of the 
ground state. To examine the influence of noise when seeking an optimal control field, first 
consider the case of aiming for a high yield (i.e., perfect control with the target yield set at 
Ot = 100% in Eq.(J2J|). Table II shows how the optimal control fields fight with different 
types of noise in this case. In order to reveal the separate contributions of noise upon the 
optimally controlled dynamics, Table II shows the yield from noise alone (o \E%\t)] 



N 

without the control field being present. Here En (t) is the z-th member of a pure field noise 
ensemble: 

E$> (t) = S(t) ^ cos (W + Tjf) • (10) 
I 

The control field is determined optimally in the presence of noise; as a reference, the yield 
O [Eoit)] of each optimal field without noise is presented. The second and third rows of Table 
II show that the optimal control can fight against amplitude noise by increasing the fluence 
F , although the effect of noise is always deleterious. The variance naturally increases as 
the amplitude Ta rises. The last three rows of Table II reveal that optimal control fields can 
be found that are fully robust to strong phase noise with small yield variance. This result 
may appear to be surprising, as proper adjustment of control field phase is often viewed as 
necessary for attaining a good quality control yield c 32j]. Table III further examines this issue 
where it is found that robustness at any level of noise assures robustness at other levels. 
However, this result alone does not reveal the degree of average coherence maintained in the 
dynamics. The shot-to-shot signal averaged density matrix of the quantum system can be 
defined as: 
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^) = ^X> w w ( lla ) 

i=l 

pV(t) = \i,([EW(t)],t))(iP([E®(t)],t)\ (lib) 

Here, lip ( [.E'W(i)] , t)) is the state driven by the z-th field in Eq.fjSJ) from the ensemble. The 
density matrix is further decomposed into its coherent and diagonal portions, R c (t), Rd(t), 
respectively, 



Rdt) = Y, M)\ 2 ; Rd{t) = Y M)\ 2 ( 12a ) 

k^j k 

R(t) = Y \PkM 2 = Rc(t) + R d (t) = Tr [p 2 (t)} . (12b) 

k,j 

Figure 3 plots these measures of the density matrix and the state populations of the system 
driven by the field in the last row of Table II with strong phase noise Tg = it. Figure 
3(al) shows that the high level of phase noise totally destroys the coherence of the system 
(R c (t) ~ 0) while Figure 3(bl) shows that large coherence is built up by the same control if no 
noise is present. The state populations are also shown in Figures 3(a2) and 3(b2) indicating 
in both cases that the target goal is reached via a simple four-step ladder climbing pathway 
which is evidently not sensitive to phase. The result in Figure 3 needs to be understood 
in the context that p^(t), for the arbitrary i — th member of the ensemble, will generally 
have a significant degree of coherence along its evolution. The details of the shot-to-shot 
coherence changes rather randomly such that R c (t) ~ 0, but the noise average (Trp 2 ) could 
contain a coherent contribution. Further insight into this situation will be analyzed upon 
examination of model 2 in Section 4.2. 

If we accept a low control yield outcome, very different control behavior is found in the 
presence of noise. The results from optimizing Eq.fjHJ) are shown in Table IV with various 
levels of amplitude and phase noise. In all cases the target yield is Ot = 2.25%. Table IV 
shows that low field amplitude noise (r^ < 0.01) has little impact on the control outcome. 
Additional amplitude noise becomes more helpful indicated by the reduced control field 
fluence. In the domain 0.06 < V a < 0.08, the control outcome is much larger than the sum 
of the yield from the control field and noise alone. For example at Ta = 0.07 the noise and 



S 



optimal field alone, respectively, produce yields of 0.26% and 0.16%. But, the same field 
operating in the presence of the noise produces a yield 2.11%. This behavior indicates that 
the field is cooperating with the presence of noise to more effectively achieve the posed goal. 
The left and right variance, and cr^, respectively, indicate a very asymmetric distribution 
having a long tail on the high population side. Again, we also find several solutions that are 
very stable to phase noise. The collective results for amplitude noise are plotted in Figure 
4. 

Figure 5 shows the evolution of the density matrix measures R c , Rd and R in Eq.([T2"|) at 
different field noise amplitude Ta levels. With strong noise there is little surviving coherence 
indicated by R c , except in the time around t~ 100 fs. In order to further explore the robust 
nature of the controlled dynamics in this low yield regime, Table V presents the outcome 
of applying noise at the amplitude T' A to a field optimally determined at noise amplitude 
T^. The yield for E^P with T' A = Ta is the outcome for the field optimally determined with 
the noise at Ta- It is evident that once a field cooperates with noise at level Ta, then it 
cooperates with noise at any other level. We may conclude that at low yield, cooperation 
can be arranged, but at the expense of coherence. 

B. Model 2 

Model 2 extends the Hamiltonian of model 1 in Figure 1(a) to allow for two quanta 
transitions: uj 02 = 2.692, U13 = 1.942, U24 = 1.314 in rad-fs -1 , with transition dipole 
elements: /i 02 = -0.1079, /i 13 = -0.1823, fi 24 = -0.2786, in 1.0 x 10~ 30 C-m. Table VI 
shows that for a high target yield of Ot = 100%, the control field can again effectively fight 
against amplitude noise by increasing the fluence F Q . But, as with the results in Table II, 
increasing the amplitude of the noise makes this battle harder with the final average yield 
decreasing and the variance increasing. The case of pure phase noise at Tq = 1.0 shows 
sensitivity to its presence with the yield falling to 68%. The last row of Table VI also shows 
a solution that is extremely robust to phase noise with reduced variance at a surprisingly 
small fluence(i.e., the fluence is even smaller than for the case Ta = 0.0.) This latter behavior 
is likely the result of there being many solutions to the optimal control search process with 
the algorithm finding alternative solutions under different circumstances. 

Table VII further analyzes the robustness of the last two phase noise contaminated fields 



9 



in Table VI. When the field Eq P (i.e., determined with no noise) is applied with a high 
level of phase noise, at T e = 2.0, there is poor robustness with the final yield dropping to 
24.2%. However, for the yield optimally determined at the high noise level T' e = 2.0, there is 
robustness to lower levels of noise. The power spectra for the field optimized at Tg = 0.0 and 
2.0 are shown in Figure 6. Comparing the two plots we see that E^^i) drives the system 
along a path of single quantum transitions while E^ {t) drives the system as well with the 
additional two quanta transition 2 — > 4. The presence of the extra two quantum transition 
can set up a delicate situation which is sensitive to phase noise, as Table VII indicates. The 
robustness of the single quantum path for Tg = 2.0 is analogous to what was found in model 
1. 

Two tables similar as Table IV and V but for model 2 can be constructed to show strong 
cooperation between optimal control field and noise when target yield is low (Ot = 3.65%), 
and the data can be plotted to form a figure analogous to that of Figure 4. 

C. Model 3 

This model is same as model 1, but amplitude noise is only present in the 1 transition. 
Thus, noise alone cannot drive any population into the target state. There is also no phase 
noise in any of the transitions. Table VIII shows the results from different levels of <-> 1 
amplitude noise for a low yield target of Ot = 10%. There is little evidence of cooperation 
with the noise, except in the last row with T A = 0.2. In this case, the field fluence is low and 
the yield from the field alone is almost zero. But, the desired target is achieved when the 
field cooperates with the noise. Figure 7 shows the power spectrum of the optimal field for 
this case and it is evident that the 0^1 transition is absent. The mechanism of cooperation 
is clear: noise drives the system from level to level 1, then the control field drives it to 
level 4. The remaining cases in Table VIII do not exclude the possibility that fields could 
be found that also cooperate with other levels of noise. The result at T A = 0.20 shows that 
such noise cooperative solutions exist. 
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D. Model 4 

The fourth model is the more complex two-path system in Figure 1(b). In this model, 
population can be transferred to the target state |4) along two separate pathways. The 
transition frequencies and dipoles of the left path are same as that of model 1; the right 
path has the distinct transition frequencies: uuo'v — 2.153, ojyy = 1.346, U2'3' = 0.345, 
co> 3 ' 4 / = 0.162 in rad-fs -1 and dipole elements: jUp'r — 0.6526, fMy^ = 0.7848, /i2'3' = 0.9023, 
A*3'4' = 1.0322 in 1.0 x 10~ 30 C-m. For simplicity, along both paths only single quanta 
transitions are allowed. 

The results in Table IX show cooperation between the field and the noise in model 4 for 
a low target yield of Ot = 2.25%. Figure 8 shows the power spectra of the three optimal 
fields in Table IX. Panel (a) indicates that in the case of no noise the field primarily drives 
system along the right path. When noise is present in the left path (panel (b)), the optimal 
field drives the system along the left path in order to cooperate with the noise. Similarly, 
when the noise is in the right path (panel (c)), the optimal field drives the system along 
the right path. An interesting feature of case (c) is the lack of intensity for the last step 
1 3') — > 1 4') along the path. Evidently the optimal field fully draws on the noise for this step 
as a special case of cooperation. The presence of the fluence term in the cost function favors 
such behavior provided that the target goal can be met. 

These results clearly show that efficiently achieving the yield is the guiding principal 
dictating the nature of the control field and mechanism in this case. This behavior reflects 
the nature of the cost function in Eq.(J7J), which places no demand on the coherence of the 
dynamics, but only asks that the control be efficiently achieved with a > 0. If coherence 
maintenance was also included in the cost function the field would be guided accordingly, in 
keeping with the ability to meet the newly posed objective. 

V. CONCLUSIONS 

The impact of amplitude and phase noise upon quantum control is explored in this 
paper with particular emphasis on these being strong. Numerical simulations from several 
cases suggest that control fields can be found that either cooperate with or fight against 
noise, depending on the circumstances. In the case of low target yields, the control field 
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can effectively cooperate with amplitude noise to drive the dynamics while minimizing the 
control ffuence. In this regime pure phase noise has little impact, if the pathway to the target 
involved single quanta transitions. But, phase noise was found to be disruptive if there are 
multiple interfering pathways. Optimal solutions were also found that exhibited stability 
with respect to phase noise while maintaining high yields. In general, when successfully 
either fighting or cooperating with significant noise, the expense paid is a measured loss of 
dynamical coherence. In some applications a lack of full coherence can be tolerated. 

The interaction between the noise and field driven dynamics is generally a highly non- 
linear process. The noise can be constructive 



quantum dynamics, improve 



24 



23| or destructive |9] in the manipulation of 
25| or reduce [38( the convergence rate of searching pro- 



cesses, and this disparate behavior makes it difficult to precisely identify the impact of noise 
under various circumstances. This work explored several special cases to reveal what may 
happen when searching for optimal controls in the presence of field noise; a mathematical 
analysis of the observed behavior may also be carried out (3]. The ability to find robust, 
and even cooperative solutions, is most encouraging. These results are consistent with a 
recent analysis [7] arguing that optimally controlled quantum dynamics has a special inher- 
ent degree of robustness. Although, each case will likely have its own unique characteristic 
behavior, the findings are encouraging for the ability to manipulate quantum dynamics even 
with severe noise being present. It would also be interesting to explore these issues in the 
laboratory by selectively introducing measured levels of noise into the control field. 
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Table I. Population distribution of the single path model 1 when the system is only 
driven by amplitude noise. 









Population 


in the state 




■p a 
J- A 





1 


2 


3 


4 


.10 


41.6% 


32.1% 


17.9% 


6.1% 


2.25% 


.09 


44.0% 


34.4% 


16.0% 


4.4% 


1.24% 


.08 


47.8% 


35.9% 


12.9% 


2.8% 


0.63% 


.07 


54.1% 


34.1% 


9.8% 


1.7% 


0.26% 


.05 


69.8% 


26.0% 


3.9% 


0.35% 


0.026% 


.03 


86.8% 


12.5% 


0.68% 


0.023% 


0.0006% 


.01 


98.4% 


1.6% 


0.0096% 


3.6xl0~ 5 % 


1.0xl0~ 7 % 


.00 


100% 


0% 


0% 


0% 


0% 


a 


The nominal field is zero: A® = 


and there ] 


s no phase noise: 7g ; 


= 0.0; T A is the 


maximum noise amplitude: 


it < ?A, 


in Eq.flU). 
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Table II. Optimal control field fighting against amplitude and phase noise with model 1 
for the high target yield Ot = 100% 



r^r^ a 


<° 




/ N 


F C 


O [E (t)) d 


F e 

L n 




a N 


n n 

U,U 




0.00% 




u.uuu 


Q8 7% 
yo. i /o 


n nn 

U.UU 


yo. 1 /o 


U /0 


0.05,0 




0.026% 


1.73 


98.6% 


3.3 x 10~ 3 


89.4% 


5.31% 


0.10,0 




2.25% 




12.3 


95.8% 


1.3 x 10~ 2 


73.6% 


13.1% 


0,1 








0.065 


97.8% 




97.1% 


0.65% 


0,2 








0.065 


95.9% 




96.2% 


1.04% 


0,7T 








0.065 


97.7% 




96.0% 


1.08% 



a Fa, r# are the maximum of the amplitude and phase noise, respectively: 



(0 



< r, 



Te. 



£# J (t) ^ : Yield from alone (Eq. ()10|l) without a control field. 
c Fluence of the control field (Eq.((3b[)). 

d O [i?o(t)]: Yield arising from the control field (Eq.tj2ajl) without noise. 
e Fluence of the noise. 

f (O \E^'{t)\) N '- Yield from the optimal control field in the presence of noise (Eq.(JHJ)] 
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Table III. Yield from optimal control fields in the case of model 1 with different levels of 
phase noise for the high target yield of Or = 100%. 









a 




V 
1 e 


E op 




■pop 
E 2 .o 


E o P 


0.0 


98.7% 


97.8% 


95.9% 


97.7% 


1.0 


97.2% 


97.1% 


96.4% 


96.7% 


2.0 


95.4% 


96.1% 


96.2% 


96.1% 


7T 


95.1% 


95.9% 


96.1% 


96.0% 



a Yield at the field E° P e (t) with phase noise: 7^ < r' e . Here E° p g (t) is the field for which 
the yield was optimized (i.e., the diagonal elements in the Table) at the phase noise level 

rv 
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Table IV. Optimal control with both amplitude and phase noise for model 1 with a low 
objective yield of Q T = 2.25%. 



r A ,r e a 






> 

/ N 


F C 




O [E (t)] d 


F e 






a N (a N ,a^) 9 


0.10,0 




2.25% 




0.00 




0.00% 


1.33x10" 


-2 


2.25% 


4.7%(2.0%,10%) 


0.09,0 




1.24% 




3.1x10" 


3 


2xl0" 6 


1.02x10" 


-2 


2.09% 


5.8%(1.8%,12%) 


0.08,0 




0.63% 




6.1x10" 


3 


2.5xl0" 5 


8.53x10" 


-3 


2.17% 


6.5%(2.0%,15%) 


0.07,0 




0.26% 




8.4x10" 


3 


0.16% 


6.53x10" 


-3 


2.11% 


5.9%(2.0%,13%) 


0.05,0 




0.026% 




9.9x10" 


3 


1.34% 


3.33x10" 


-3 


2.08% 


4.7%(1.8%,9.3%) 




6.0xl0" 4 


% 


X . \J\J A 1U 


-2 


1 87% 


1 90x10" 


-3 




9 6%f1 4% 4 1 %) 


0.01,0 


1.0xl0~ 7 % 


1.11x10 


-2 


2.12% 


1.33x10" 


-4 


2.17% 


0.8%(0.7%,1.0%) 


0,0.5 








1.09x10 


-2 


2.18% 






2.17% 


0.01% 


0,1.0 








1.10x10- 


-2 


2.21% 






2.18% 


0.03% 


0,2.0 








1.10x10 


-2 


2.19% 






2.17% 


0.03% 


0,0 








1.09x10 


-2 


2.20% 






2.20% 


0.00% 



a ,b,c,d,e,f re f er to Table II. 



9 CTjy, a^j are left and right standard deviation, respectively, around the mean in 
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Table V. Yield from optimal control fields with different levels of amplitude noise for 
model 1. 









tp[ 




a 








r" 


E op 


-^0.09 


E op 


E op 


E op 

-^0.05 


E op 

-^0.03 


E op 


E op 
-^o.oo 


0.10 


2.25 


2.98 


4.15 


4.35 


4.64 


4.74 


5.25 


5.50 


0.09 


1.24 


2.09 


2.96 


3.80 


3.99 


4.11 


4.79 


4.53 


0.08 


0.63 


1.27 


2.17 


2.92 


3.28 


3.60 


4.16 


3.76 


0.07 


0.26 


0.71 


1.48 


2.11 


2.90 


3.18 


3.48 


3.46 


0.05 


0.026 


0.17 


0.55 


1.1 


2.08 


2.50 


2.91 


3.36 


0.03 


6.0xl0~ 4 


0.023 


0.14 


0.44 


1.58 


2.06 


2.38 


2.45 


0.01 


1.0xl0~ 7 


0.0014 


0.016 


0.19 


1.38 


1.88 


2.17 


2.21 


0.00 


0.00 


0.00020 


0.0025 


0.16 


1.34 


1.87 


2.12 


2.20 



a Yield at E° p (t) with noise bounded by It^J < T' A and for 7^ = 0. Here E^ P A (t) is 
field for which the yield (i.e., the diagonal elements in the table) was optimized with noise 
T' A = T A . The yields are expressed in percent. 



19 



Table VI. The yield from optimal control fields fighting with amplitude or phase noise in 
model 2 with the high yield target objective of Ot = 100%. 



r^4,r# a 


\ 




In 


O[E (t)\ d 


F n e 


(O L f 

\ L v / J 1 N 




0.0,0.0 




0.00% 


0.28 


99.1% 


0.00 


99.1% 


0.00% 


0.05,0.0 




0.13% 


1.53 


98.2% 


5.83xl0- 3 


87.2% 


7.4% 


0.1,0.0 




3.65% 


3.27 


93.7% 


2.33xl0- 2 


66.5% 


14.5% 


0.0,1.0 






0.33 


80.4% 




68.0% 


13.1% 


0.0,2.0 






0.065 


96.5% 




91.6% 


4.2% 



,b,c,d,e,f re f er to Table II 
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Table VII. Control yields with optimal fields at different phase noise levels for the high 
yield target Or = 100% in case of model 2 



V 
1 e 




(o [K p M) K a 




E op 

^0.0 


■pop 


E op 


0.0 


99.1% 


80.4% 


96.5% 


1.0 


42.9% 


68.0% 


95.3% 


2.0 


24.2% 


37.5% 


91.6% 



a refer to Table III 
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Table VIII. Control outcomes from optimal fields for model 3 with a low target yield of 



T = 10%. 


T A a I Pi 
\ 




) b Fo c 

1 N 


O [E (t)} d 


(O [EM(t)]) N f 


cr N 


0.00 


0.00% 


0.0176 


9.95% 


9.95% 


0.00% 


0.01 


1.67% 


0.0198 


9.91% 


9.94% 


1.39% 


0.05 


31.7% 


0.0229 


10.66% 


9.97% 


3.24% 


0.10 


61.1% 


0.0194 


9.27% 


9.99% 


8.09% 


0.20 


46.6% 


0.0144 


0.0013% 


9.90% 


6.06% 



a Amplitude noise is only in the transitions <-> 1. 

6 Population in level 1 when the system is only driven by noise. 

c,d,f re fer to Table II. 
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Table IX. Yield attained from the optimal fields for model 4 with the low yield target of 
T = 2.25%. 



-p(L) -p(-R) a 
1 A ' 1 A 






) b F C 

1 N 


O [E (t)] d 


(O mt)]) N f 


0.00,0,00 




0.00% 


0.0093 


2.19% 


2.19% 


0.08,0.00 




0.63% 


0.0064 


0.0003% 


2.17% 


0.00,0.07 




0.50% 


0.0052 


0.0002% 


2.08% 



a F\ F ( A ': Noise in the left and right paths, respectively; 
b,c,d,f re f er to Table II. 
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FIG. 1: Two quantum multilevel systems used to investigate the impact of noise on optimally 
controlled dynamics, (a) The 5 level single path system used for simulations with models 1,2,3 in 
Sections 4.1, 4.2, and 4.3, respectively. Model 1 and 3 only allow single quanta transitions while 
model 2 allows both single and double quanta transitions, (b) The double-path system for model 
4 in Section 4.4. 
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FIG. 2: Optimal control fields and their power spectra for model 1 when the field has no noise. 
The fields are found using the cost function in Ea. paf) . Plots (al) and (a2) correspond to the 
high target yield of Ot = 100% while plots (bl) and (b2) correspond to the low target yield of 
O t = 2.25%. 
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FIG. 3: The time dependence of model 1 dynamics driven by the optimal control field in the last 
row of Table II. The density matrix measures Rd(t), R c (t), R(t) in plots (al) and (bl) are defined 
in Ea.(|12[). Plots (al) and (a2) show the dynamics when the system is driven by a control field 
with noise while plots (bl) and (b2) show the dynamics of the system driven by the same field but 
without noise. The associated state populations are shown in plots (a2) and (b2). 
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FIG. 4: The control yield under various noise conditions for model 1 with the low yield target of 
Ot = 2.25% in Table V. The labeled curves show the yield for noise alone, the yield for the optimal 
field alone and the yield for the field in the presence of noise. There is notable cooperation between 
the noise and the field especially over the amplitude noise range 0.06 < Ta < 0.08. 
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0.5 
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Time(fs) Time(fs) Time(fs) 



FIG. 5: Evolution of the density matrix measures in Eq.(12) with different levels of field amplitude 
noise for the low target yield case of Figure 4. The coherence is successively destroyed by stronger 
noise, but it still survives around t~ 100 fs. 
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FIG. 6: Power spectra of the control fields for model 2 aiming at a high yield with Ot = 100%. 
Plot (a) is for the field E^ Q robust to phase noise; plot(b) is for the field Eq P which is not robust 
to noise. The robust control field in (a) avoids the extra interference arising from the two quanta 
transition 2^4. 
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FIG. 7: The optimal field power spectra of model 3. Here noise is only in the amplitude of the 
<-> 1 transition with Ta = 0.20. The optimal field cooperates with the noise by eliminating the 
ujqi frequency and exclusively relying on the noise to drive that step. The arrow indicates where 
that frequency lies. 
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FIG. 8: Power spectra of the control fields for model 4 (see panel (b) of Figure 1). The frequencies 
along the right path in Figure 1 are denoted with a prime to distinguish them from the left path, 
(a) No noise, (b) noise in the left path, (c) noise in the right path. In these low target yield cases 
the optimal field cooperates with the noise and follows the associated path. 
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